引言
- 定义:特征值全是正数的实对称矩阵为正定矩阵(positive definite matrix)。
- 类似的,若实对称矩阵的特征值均非负,则为半正定矩阵(positive semidefinite matrix)。
可能用到的概念
实对称矩阵A正定的充要条件
注意,这里的所有进行判别的矩阵都是实对称矩阵。
以下条件都是判别实对称矩阵是否正定的充要条件。
- A的所有特征值λi均为正。
- xTAx>0对所有非零向量x都成立。
- A的所有顺序主子式都是正的。
- A的所有主元(无行交换)都是正的。
- 存在列满秩矩阵R,使得A=RTR.
- A的所有主子式都是正的。
证明
(1) => (2)
对实对称矩阵A,存在正交阵Q,使得A=QΛQR。
因此,对任意非零向量x:
xTAx=xTQΛQTx=yTΛy=λ1y12+...+λnyn2>0
其中y=QTx=(y1,...,yn)≠0
(2) => (1)
因为A为实对称矩阵,因此一定可以对角化A=QλQT
∵xTAx>0∴xTQΛQTx>0
令y=QTx,因此yΛyT>0⇒∑λiyi2=0
∵∀y≠0s.t∑λiyi2=0
λi=0
(2) => (3)
由(2) => (1) => det(A)=λ1...λn>0
xTAx=(xkT0)⎝⎛Ak∗∗∗⎠⎞⎝⎛xk0⎠⎞=xkTAxk>0
∴det(xkT)det(Ak)det(kk)=det(Ak)>0
(3) => (4)
顺序主子式与主元有直接关系,第k个主元:
dk=detAk(k−1)det(Ak)>0
其中Ak是第k个顺序主子矩阵(the k-th leading principal submatrix)。
(4) => (2)
由对称矩阵的Gauss消元法得LDU分解:A=LDLT,其中对角阵的对角元为A的主元:
D=diag(d1,...,dn)
则对任意非零向量:
xTAx=xTLDLTx=yTDy=d1y12+...+dnyn2>0
(2) => (5)
A=LDLT=LDDLT=(DLT)T(DLT)=RTR
(5) => (2)
设A=RTR,则对任意非零向量x:
xTAx=xTRTRx=(Rx)T(Rx)=∣∣Rx∣∣2>0
(6) => (2)
(6) => (3) => (2)
(2) => (6)
对k阶主子矩阵Ai1,...,ik,任取x=(x1,...,xn)≠0,使其除xi1,...,xik的其余分量全为0,则
xTAx=(xi1,...,xik)Ai1,...ik⎝⎜⎜⎜⎜⎛xi1...xik⎠⎟⎟⎟⎟⎞>0
∴det(Ai1,...ik)>0
判断
因此,我们可以用以上充要条件来判断矩阵是否正定。常用的判断方法有:
- 看顺序主子式是否都大于0
- Gauss消元后主元是否都大于0
- 看特征值是否都大于0
- 找任意一个向量计算xTAx是否大于0
- 找A=LDLT分解,看R=DLT是否满秩
正定矩阵的常见性质
正定矩阵一定可逆。(行列式大于0)
设A,B为正定矩阵,则A+B也为正定矩阵。
设A为正定矩阵,则存在矩阵C,使得A=C2。
设A为正定矩阵,则矩阵A2和A−1也正定。
设A为正定矩阵,矩阵C可逆,则B=CTAC也正定。
半正定矩阵
充要条件
- A的所有特征值λi均非负。
- xTAx≥0对所有向量X成立。
- 存在矩阵R,使得A=RTR(R可能不是可逆阵)。
- A的所有主子式均非负。(注意,不是所有顺序主子式)